TECHNICAL 

MEMORANDUM 


AN ANALYTICAL SOLUTION TO PATCHED 
CONK: TRAJECTORIES SATISFYING INITIAL 
AND FINAL BOUNDARY CONDITIONS 





Bellcomm. Inc. 

955 I'ENfANT PLAZA NOWTH, S.W., WASHINGTON, D.C. 20024 


COVER SHEET FOR TECHNICAL MEMORANDUM 

title- An Analytical Solution to Patched tm-70-2011-1 

Conic Trajectories Satisfying 

Initial and Final Boundary Conditions date- November 30, 1970 

FILING CASE NOISI-J^Q 

autmor<*>- K. M. Carlson 


filing subject(s) Patched-Conic Trajectories 

(ASSIGNED *V AUTHOR.*, ,- Interplanetary Trajectorie s 

Lunar Trajectories 
Trajectory Analysis 
Celestial Mechanics 


ABSTRACT 


Analytical expressions for patched-conic trajectories 
are derived, thus providing the analyst with a tool for the 
calculation of spacecraft trajectories between two bodies and 
for understanding the behavior of these trajectories as the mis- 
sion constraints are varied. Two sets of equations are developed 
which allow an analytical solution to spacecraft trajectories 
between two large central bodies and which satisfy boundary con- 
ditions at both bodies. The derivations of these equation sets 
are based on patched-conic analysis. The independent parameters 
required are the spacecraft angular momentum and inclination 
with respect to the more massive body, and the spacecraft angular 
momentum, energy, and inclination with respect to the less mas- 
sive body, or equivalent quantities. These equation sets provide 
a good approximation to the sphere of influence trajectory patch 
point over essentially the entire range of possible trajectories. 
Having obtained the patch point, the spacecraft state vector is 
easily obtained. 
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The study of spacecraft trajectories between two 
bodies does not yield to mathematical analysis as well as 
trajectories about a single body. Consequently, the under- 
standing of such trajectories has been difficult and, to 
an extent, intuitive. This memorandum provides a major 
simplification in the analysis of such trajectories. 

The analysis contained here is based on the patched- 
conic approximation to trajectories between two large central 
bodies. In Reference 1, the author carried out a detailed 
analysis of the basic elements of the trajectory patch point 
locus on the sphere of influence. These elements were iden- 
tified and described on a geometrical basis. Using the 
geometrical insights gained in Reference 1, equations are 
derived in this memorandum which predict the location of the 
patch point using trajectory parameters from both inside and 
outside the sphere of influence. Once the patch point is 
located, the trajectory state vector can be easily calculated, 
and hence the spacecraft position may be found for any time. 
The trajectory parameters required by the patch point equa- 
tions can be thought of as being the elements of the trajec- 
tory. There are six such elements, just as there are six 
elements of a trajectory about a single body. Since the 
parameters describe features both inside and outside the 
sphere of influence, obtaining a trajectory with desired 
characteristics is considerably simplified. 

Two equation sets are derived. One set is valid 
over essentially the entire regime of possible trajectories. 
The other equation set is valid only for trajectories at the 
high energy end of the trajectory regime; however, this set 
is considerably simpler than the set with general validity. 

The basic elements of the sphere of influence 
patch point locus will be discussed first without proof. 
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The patch point equation sets will then be derived and the 
necessary proofs supplied. Finally, the validity of the 
analysis will be illustrated numerically, and some useful 
techniques and approximations will be presented. The patched- 
conic equation sets are restated in Appendix B for easy 
reference. 


Several terms which are used frequently in the 
analysis bear some explanation. The two large central bodies 
are referred to as the major body and the minor body , the 
major body being the more massive of the two. T he sphere of 
influence is a spherical region surrounding the minor body 
within which only the minor body's gravitational field is 
assumed to operate. All the rest of space is assumed to be 
influenced only by the major body's gravitational field. 
Trajectories are named descriptively. Thus, a rectilinear 
trajectory is one which rises and/or descends vertically with 
respect to its central body, i.e., the spacecraft position and 
velocity vectors are colinear. A dual rectilinear trajectory 
passes between the major and minor bodies and consists of 
rectilinear trajectory segments centered at each body. A 
single rectilinear trajectory again passes between the two 
bodies but has only one body-centered segment which is rec- 
tilinear. A non-rectilinear trajectory is one such that 
neither of its body-centered portions is rectilinear. A re- 
statement of these definitions, as well as a list of symbols 
used, is contained in Appendix A for easy reference. 

2.0 STRUCTURE OF THE SPHERE OF INFLUENCE PATCH POINT LOCUS 


The sphere of influence patch point locus for non- 
rectilinear trajectories can be constructed from the dual and 
single rectilinear trajectory patch point lo^i with the same 
energy. Figure 1 shows the dual and single rectilinear tra- 
jectory patch point loci for trajectories with identical minor 
body referenced energies. 

It is seen from Figure 1 that the dual rectilinear 
patch point locus (C) is a single point lying on the intersec- 
tion of the sphere of influence and the minor body orbital 
plane. The major-body-centered and minor-body-centered 
single-rectilinear trajectory loci (A and B) are closed curves 
containing the dual rectilinear patch point. Each point of the 
single rectilinear trajectory loci (e.g., A and B) represents 
a different trajectory plane orientation. The sizes of the 
single rectilinear loci are dependent on the angular momenta 
of the non-rectilinear trajectory segments. Angular momentum 
is, in turn, relaced to pericenter radius. Note that the 
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trajectory planes for single rectilinear trajectories paii, 
through the dual rectilinear patch point. That is, the node 
line between the spacecraft trajectory plane and the minor 
body orbital plane is the path of a dual rectilinear trajec- 
tory with the same minor body referenced energy. 

Figure 2 shows the patch point locus for non- 
rectilinear trajectories. It will be shown that this locus 
can be constructed as follows: 

1. Calculate the dual and single rectilinear trajectory 
patch point loci (A, B, and C) for the desired 
trajectory energy and angular momenta as shown in 
Figure 1. 

2. Select the minor-body-centered single-rectilinear 
patch point (B) representing the major-body- 
centered trajectory plane desired. 

3. Transpose the major-body-centered single-rectilinear 
patch point locus (A) to A' so that the interior 
point originally coinciding with the dual rectilinear 
patch point (C) now coincides with the point B 
selected in Step 2. 

Each point of the locus created in Step 3 is a good 
approximation to the patch point for a non-rectilinear trajec- 
tory with a given minor-body- centered plane orientation, the 
major body plane orientation used in Step 2, and the angular 
momenta and energy used in Step 1. In other words, the dis- 
placement of the patch point caused by going from a dual recti- 
linear trajectory to a major-body-centered single-rectilinear 
trajectory is essentially identical, in magnitude and direction, 
to the patch point displacement caused by going from a minor- 
body-centered single-rectilinear trajectory to a non-rectilinear 
trajectory. Note that this also implies that the minor-body- 
centered plane passes through the minor-body- centered single- 
rectilinear patch point. 

3*0 SPHERE OF INFLUENCE PATCH POINT EQUATION DERIVATION 
3.1 Coordinate Systems 

In the following analyses, two coordinate systems 
are used, one centered at each body. The orbital plane of the 
minor body is taken as the reference (xy) plane of both systems. 
The line joining the major and minor bodies at the time the 
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FIGURE 2 - STRUCTURE OF THE NON-RECTI LINEAR PATCH POINT POINT LOCUS 
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spacecraft crosses the sphere of influence is taken as the x 
axis for both coordinate systems. The positive x direction is 
taken as toward the body not at the system center. Thus, the 
major-body-centered positive x axis points toward the minor 
body and, conversely, the minor-body-centered positive x axis 
points toward the major body. For both coordinate systems, 
the positive z axis points in the same direction as the angular 
momentum vector of the minor body in its orbit around the major 
body. The positive y axis completes the right-handed triad. 

The transformation relationships between these systems are: 


X 1 = 


-x 


.1 


*1 = - y 2 > 


( 1 ) 


Z 1 " 


+z 


V 


3 . 2 Dual Rectilinear Trajectory Patch Point 

Figure 3 shows the trajectory geometry required for 
dual rectilinear trajectories. From the geometry of the figure, 
the minor-body-centered state vector at the sphere of influence 
is 


R i ■ R i 


COSX DRT eOSfi DRT 


sinx DRT eo8 6 DRT 


sine 


DRT 


and, since R^ and V, are colinear. 


V 1 - 1 * 


C08 X drt C0>6 drt 


SinX DRT COS8 DRT 


sing 


DRT 


( 2 ) 


(3) 




dual RECTILINEAR 
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where k determines the direction of and obeys the following 
rules: 


k * +1 for trajectories from the minor body to the 
major body. 


k = -1 for trajectories from the major body to the 
minor body. 


The major-body- centered state vector at the sphere of influence 
is obtained from the vector sums 


R 2 = 


R 1 + R o 


and 


( 4 ) 


V 1 * 


( 5 ) 


4 ^ 

where and must be transformed to the major body coordinate 

system using Eqn. (1) . The state vector of the minor body 
referenced to the major body is given by 


and 


1 
0 
0 . 

-cos * 0 

+8in 4 ft 


( 6 ) 


( 7 ) 
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In order to obtain a major-body- centered rectilinear trajec- 
tory , the angular momentum with respect to the major body must 
be equal to zero, that is 


R 2 x V 2 = 0. (8) 

Writing Eqn. (8) in terms of Eqns. (2) through (7) , 


- R l V o sin6 DRT 8ir % 

- < R l V o C ° 8 *o + Vl k > 8in6 DRT _ Q (9) 

R o V o sin *o " R lV in *o C08X DRT C08B DRT 

- (R^cos^ + R^k) sinx DRT cos8 DRT 

Setting the first and second components of the angular 
momentum equal to zero yields* 

sinBjjRT = 0. (10) 


That is, dual rectilinear trajectories must lie in the orbital 
plane of the minor body. 

Setting the third component of angular momentum equal 
to zero (and setting cos& DRT * 1) 

( R l v o cos *0 + R o V l k)8inX DRT + R l V 0 8in ^ 0 CO8X DRT “ R o V o 8in ^0 * °* 


*An alternative solution is 
sint c - 0 

* nd R l v o 

V l k * K — co **o 

However, this solution is of academic interest only, since 
it corresponds to a rectilinear trajectory of the minor 
body with respect to the major body. 
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Equation (11) has the form 


Asino + Bcosct + C » 0 


which can easily be solved, yielding 


sina = 


-AC±B YA 2 + B 2 - C 2 


A + B 


We now have the equations of the dual rectilinear 
trajectory patch point: 


sin6 DRT * 0 


sinx. 


-AC*B VA 2 + B 2 - C 2 


A + B 


where 


A = R.V cos* A + R V,k 
1 O O O 1 

B = R.V sin* 
l o o 


C «* -R v _sin<j>_ . 
o o r o 


Thus, e DRT is constant (0°) and X DRT is a function of Vjk only. 

Due to the sign ambiguity on the radical, the * DRT 
equation yields two solutions except when 


C 2 * A 2 + B 2 
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It can be shown from purely geometrical considerations that 
Eqn. (15) is the condition for a major-body- centered trajec- 
tory which is tangent to the sphere of influence. Since it 


2 2 2 

is not possible to have C < A + B without introducing 
imaginary values in the x DR(T , equation, Eqn. (15) also imposes 

a lower limit on the value of V^. Solving Eqn. (15) for V^, 


V 


1 


R. V COS A 
1 o T o 


+ ->/R 


£ 



V in *o 


(16) 


Equation (16) represents the minimum energy for a dual recti- 
linear trajectory, and 


sinx 


DRT 



(17) 


is the patch point for the minimum energy trajectory. For 
trajectories with energies greater than that given by Eqn. (16) , 
two values of x DRT exist, representing trajectories piercing 

the sphere of influence on either side of the tangent trajectory. 
Thus, the trajectories are grouped by 

X DRT < sin (" C ) 

or (18) 

(b) X DRT > sin_1 (" I ) 

The sign operating on the radical is minus for case (a) and 
plus for case (b) . x DRT lies in the first or second quadrant 

for k = +1 and in the third or fourth quadrant for k * -1. 

3 . 3 Major- Body-Centered Single-Rectilinear Trajectory Patch 
Point Locus 


Figure 4 shows the trajectory geometry required for 
majcr-body-centered single-rectilinear trajectories. From the 
geometry of the figure, the minor-body-centered state vector 
at the sphere of influence is 


t 





SPHERE OF. 
influence 
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C08Q SRT1 C08y SRT 1 + sinfi SRTl COSi l 8in Y S RTl 

Sinft SRTl COSY SRTl “ COSft SRTl COS1 l SinY SRTl 
siniiSiny SRT1 


(19) 


(-cos^j^) (cosf2 SRT1 cosy SRT1 + sin£i SRTl cosi l sinY SRTl * 

+ (sin^) (cosfl SRT1 siny SRT1 - sinflg^^osijCosyg^) 


V, 


(-cos^) (sinn SRTl cosy SRT1 - cosng RTl co8i lS iny SRT1 ) 

+ (sin^) (sim2 SRT1 sinyg RT1 + cos£ isRTl cosi, l co8Y SRTl^ 


( 20 ) 


- cos^ 1 sini 1 siny SRT1 - sin^sinijcosyg^ 


At this point the definitions used for inclination and 
node line require clarification. Three things are required to 
define the trajectory plane orientation: the angle between the 

trajectory plane and the reference plane (inclination) , the loca- 
tion of a node, and whether this node is the ascending or de- 
scending node. The standard procedure is to select the ascending 
node and measure inclination as positive. However, it is conven- 
ient here to use the node line nearest the pierce point and assign 
an algebraic sign to the inclination to show whether this node is 
ascending or descending. The convention chosen is that a positive 
inclination should produce a positive patch point latitude for 
trajectories from the minor body to the major body. This conven- 
tion requires ^ to be positive if the nearest node is descending, 

and negative if it is ascending. The symbol ft is used to denote 
the longitude of the nearest node rather than its conventional use 
denoting the longitude of the ascending node. 

Returning to the analysis, in order to obtain a major- 
body-centered rectilinear trajectory, the angular momentum must 
equal zero, that is 


-► 



X V, 


0 . 


( 21 ) 



t 
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Rewriting Eqn. (21) in terms of Eqns. (4) through (7) and 
Eqns. (19) and (20) yields 


R 1 V 1 sinfi SRT1 sini 1 sin4> 1 - R 1 V 0 sini 1 siny SRT1 sin^ o 
- R l v 1 cos0 SRT1 sini 1 sin4 1 - RiV ini l sinY SRTl cos *o 

+ R q V 1 (sinij^sinyg^j^cosi}.;,^ + sini 1 coSY SRT1 sin$ 1 

R l V l cos *l S * n< *'l 

+ R^Ksin^) (-C0sfl SRT1 C08Y SRTl - sinfi SRTl COSi l SinY SRTl ) 
+ (cos* Q ) (-sinO SRT1 cosY SRT1 + cosn S RTl cosi l sinY SRTl ) J 
+ r^Kcos^) (+sin« SRTl cosY SRTl - cosn S RTl cosi l sinY SRTl ) 

- (sin^j^) (sinfi SRTl sinY SRT1 + cosOg^cosi^osYg^)] 

= 0 

Setting the first component of the angular momentum 
equal to zero 


V^sin$^ 

sir "'sRTl " V 0 silH" sln0 SRTl' (23: 

which gives the patch point locus displacement from ^g R ^ i n 
terms of n SRT1 - Note that y srt1 is not a function of 

Now, assuming the miner body pericenter radius is 
small compared to the sphere of influence radius, then sint^ 

has a small value. If the minor body orbits the major body in 
a low eccentricity orbit, sin$ Q will be approximately one. 


( 22 ) 
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Then, since V Q and are typically of the same order of 
magnitude, and since 4» ^ grows smaller as grows larger, 
Eqn. (23) shows y srt1 to be a small angle. We may then 
make the approximations 


cos^ % -k 


and 


(24) 


cosy srt1 £ 1. (25) 

Using these approximations and Eqn. (23) , and setting 
the second component of angular momentum equal to zero? yields 


(R l V o COS *o + R o V l k) sinn SRTl + R l V o sin *o COsn SRTl - R o V o sin *o = °‘ 


Equation (26) is identical to the dual rectilinear patch-point 
equation, Eqn. (11). Consequently, 


^SRTl % X DRT* (27) 


Thus, it has been shown that the node line of a major-body- 
centered single-rectilinear trajectory is closely approximated 
by the path of the associated dual rectilinear trajectory. 

Note that this also implies that major-body-centered single- 
r ectilinear trajectories do not exist for energies below the 
minimum dual rectilinear trajectory energy of Eqn. (16) . 

From basic spherical trigonometry relations. 


and 


sxn6 SRTl " sini l siny SRTl 


tanA XgRTi = cosi 1 tan YsRT1 


(28) 


(29) 


(26) 
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where -* SRT1 is the longitudinal distance from fi SRT ^ to the 
single rectilinear trajectory patch point. Note that while 

I 

i-^ is the trajectory plane inclination for major-body-centered 

single-rectilinear trajectories (and hence identical to i-^) , 

this is no longer true when the locus is transposed to form the 
non-rectilinear locus, and a more general definition is required. 

I 

The general definition of i^ is the angle between the minor-body- 

centered trajectory plane and the plane containing the radius 
vector to the center of the Yg RT ^ locus and whose line of nodes 

with the reference plane is perpendicular to this radius vector. 

Summarizing the equations for the major-body-centered 
single-rectilinear patch point, 


and 


where 


and 


sinSgRTi = sini’sinY SRT1 " 
X SRT1 = X DRT + AX SRT1 


I 


tanAXg RT1 = cosi^tany 


V^sin^ 

siny RRTl = V sin* 
o o 


SRT1 


sinX 


DRT 


(30) 


and 'DRT is found from E< I n * d4). Note that i 1 carries a plus 

sign for trajectories whose patch point lies above the orbital 
plane of the minor body and a minus sign for trajectories whose 
patch point lies below the orbital plane of the minor body. 

3.4 Minor-Body-Centered Single-Rectilinear Trajectory Patch 
Point Locus 


Figure 5 shows the trajectory geometry required for 
minor-body-centered single-rectilinear trajectories. Note that 




MINOR-BODY-CENTERED SINGLE-RECTILINEAR TRAJECTORIES 
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the state vector is defined using the patch point coordinates, 
since there is no defined minor- body- centered plane. From the 
geometry of the figure, the minor-body-centered state vector 
at the sphere of influence is 



COSX SRT2 COS6 SRT2 

sinX SRT2 COS6 SRT2 

sin6 S RT2 


(31) 


and, since the minor-body-centered radius and velocity vectors 
are colinear, 



COSX SRT2 COS6 SRT2 

SinX SRT2 COS8 SRT2 

Sin6 SRT2 


(32) 


where k = +1 for trajectories from the minor body to the major 
body 

k = -1 for trajectories to the minor body from the major 
body. 

For this case, a non-rectilinear major-body-centered 
trajectory is required. Hence, 

R 2 x V 2 - h 2 . (33) 

The angular momentum vector may be written 
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sini 2 S in !! SRT2 

-sini 2 cos(l SRT2 

cosi 2 


(34) 


Combining Eqns. (4) through (7) and Eqns. (31) and (34) 


- R lV ine SRT 2 Sin *o 
-R 1 v 0 sinB SRT 2 cos * 0 - VV inS SR T 

- R lV cosA SRT2 COS6 SRT2 Sin *o + SinX SRT2 cos6 SRT2 c ° s *o ) 
-R 0 kV 1 sinX SRT2 c°se SRT2 + R 0 V 0 sin* 0 


= h 


2 


sini 2 sinfi SRT2 

-sini 2 cosn SRT2 

coai 2 


(35) 


From the first component of the angular momentum 

equation 


h 2 8ini 2 8intl SRT2 
8in6 SRT2 R 1 V 0 sin4 0 


(36) 


From the second component of the angular momentum equation , we 
find that 
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. h 2 3ini 2 cosf! SRT2 
sint5 SRT2 R,V cost + R kv, * 

X O O O 1 

Combining Eqns. (36) and (37) 


(37) 


R.V sini 

^ _ 1 o o 

tan “SRT2 " ‘ R.V cos* + R„kV, * 

1 O O O 1 

Note that fig RT2 i® independent of either hj or ij. 


(38) 



FIGURE 6 - GEOMETRY OF MINOR-BODY-CENTERED SINGLE-RECTILINEAR TRAJECTORY FLANE NODE 
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This being so, from Figure 6 one can write 


siR (n SRT2 + n SRT2 ) " ^ SinfJ SRT2 


( 39 ) 


provided the node line intersects the sphere of influence. 
Combining Eqns. (38) and (39) 

< R o kV l + *l V o COS *o ) sinfi SRT2 + R l V o sil % CCSfi SRT2 ’ R oV in *o = °* 


Equation (40) is identical to the X nB „ equation, E'tn. (11); 
hence, DRT 


fl SRT2 * X DRT* (41) 

Thus, it has been shown that the node line of a minor-body- 
centered single-rectilinear trajectory is the path of the asso- 
ciated dual rectilinear trajectory for energies above the minimum 
dual rectilinear trajectory energy. However, it is not implied 
in the derivation that minor-body-centered single-rectilinear 
trajectories do not exist for energies below this value. In 
fact they do exist for such energies. In these canes, the major- 
body-centered node line does not intersect the sphere of influence. 

From the third component of the angular momentum 

equation 


(R o kV l + R 1 V 0 C08 *0 ) sinX SRT2 * R l V 0 8in *0 C08X SRT2 


R^V^sin* - h~cosi 0 
o o T o 2 2 m 0 

eo,8 SRT2 

Again this is an expression of the form 

Asina + Bcosa + C - 0 


(42) 


and Eqn. (13) constitutes the solution. 


(40) 
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We now have the equations for the minor-body-centered 
single-rectilinear patch point, 


where 


h ^ s ini 2 sin & SRT2 
sin6 SRT2 = R 1 V 0 sin(f 0 

R, V sin* 

1 o lo l 

tann SRT2 ' R 1 V o cos^ o + R Q kv l > 


sinx 


AC 


± B /"a 2 + B 2 - C 2 


SRT2 


A 2 + B 2 


(43) 


A = R 0 kV x + R^cos^ 

B = R^sin^ 

R V sin<* - h-cosi- 
r _ O O T Q 2 £ 

cose SRT2 


The 


angle 


X SRT2 


lies in the first or second quadrant for k = +1 


and in the third or fourth quadrant for k = -1. The convention 
that a -ositive inclination produces a positive patch point lati- 
tude f^. .flinor body to major body trajectories requires a posi- 
tive ±2 if the nearest node is ascending and negative ij if it 


is descending. 


Because of the radical in the x gRT 2 equation, 
A 2 + B 2 - C 2 - 0 


(44) 


represents a limit on minor-body-centered single-rectilinear 
trajectories. Solving Eqn. (44) Jor 


-kR.V cosifc + 
1 o o 


V, = 


vf 


_ - h 0 cosi^ 
0 2 2 


cos6 SRT1 


-) - <v, 


o ain *o ) 


(45) 
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Equation (45) imposes a lower limit on the energy of minor- 
body-centered single- rectilinear trajectories and 


rinx 


A 

SRT2 “ C 


( 46 ) 


is the patch point for the minimum energy trajectory. 


As before, two solutions exist for ^ srt2 f° r energy 

levels above that given in Eqn. (45) . Both solutions are valid 
and are grouped by: 


(a) 


or 


(b) 


X SRT2 < Sin I cl 

► 

X SRT2 > Sin | D 


(47) 


The sign operating on the radical in the x SRT 2 equation is minus 
for case (a) and plus for case (b) . 


3.5 Non-Rectilinear Trajectory Patch Point Equations (High 
Energy Set ) 

It was stated at the beginning of the analysis that 
the non-rectilinear patch point locus can be obtained by 
properly combining the dual and single rectilinear patch point 
loci. Equations for these loci have now been derived, so the 
general locus may now be constructed. Specifically, the major- 
body-cantered single-rectilinear patch point locus is to be 
transposed so that the interior point originally coinciding 
with the dual rectilinear patch point now lies on the minor- 
body-centered single-rectilinear trajectory patch point. 

That is. 


and 


6 nrt 


X NRT 


6 SRT2 + 3 
X SRT2 + AX 


SRT1 


SRT1 




(48) 
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where the values for 0 SRT1 / 0 SRT2 ' AX SRT1' and X SRT2 are found 

I 

from Eqns. (30) and (43) . The value of i used in Eqn. (30) 

1 I 

must be obtained using the definition of i, given in Section 

I 

3.3. Spherical trigonometry easily shows that i^ may be ob- 
tained from the true inclination, i^, by the relationship 


cos 



cosi^ 

COse SRTl 


(49) 


An alternative method for determining & NRT and 

is to calculate the angular distance from the minor-body- 
centered node line to the minor-body-centered single-rectilinear 
trajectory patch point and add y srt ^ to this. The non- 
recti linear patch point coordinates may then be calculated 
using this sum and i^. The resulting equations are 



J 


While the two formulations give slightly different 
results, there does not appear to be any advantage of one over 
the other. Therefore, the slightly simpler Eqns. (48) are used 
here. 

For reference purposes, Eqns. (48) will be called the 
high energy equation set- The reason for this selection will 
become apparent. 

The independent parameters in the high energy equa- 
tion set are: 
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minor-body- centered 
inclination (i-^) 

velocity at the sphere of influence (V^) 

trajectory flight path angle 

ma jor-body-centered 

angular momentum (l^ ) 


inclination (ij). 

The five independent parameters and the time of piercing the 
sphere of influence may be considered the elements of the 
trajectory, just as there are five elements to Keplerian tra- 
jectories which, along with the time of pericenter passage, 
determine two body trajectories. 

The elements and <J> ^ may, if desired, be restated 

in terms of other parameters using the conic equations. For 
example , 



( 51 ) 


( 52 ) 


obtains and 4,^ from the semi-major axis and the radius of 
pericenter for the minor body. Means of obtaining hj are dis- 
cussed in Appendix C. The parameters R Q , V Q , and $ 0 are obtained 

from the minor body ephemeris and are calculated at the time 
the spacecraft pierces the sphere of influence. 

The results of a sample trajectory problem are 
presented to illustrate the validity of the concepts developed 
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so far. Consider a trajectory from the Moon to the Earth, such 
as a return from an Apollo mission. The selenocentric parameters 
selected are; 


periselene altitude 


60 n. mi. 


inclination 
semi-major axis 


160° 

8 x 10 6 feet 


and the geocentric parameters are: 

inclination 30° 

11 2 

angular momentum 7.7 x 10 ft /sec. 


The radius of the Moon's sphere of influence is taken as 
8 

1.8 x io feet. The Moon is assumed to be at apogee, a dis- 

Q 

tance of 13.34 x io feet from the Earth. Using standard conic 
orbit equations 


and 


= 4854.4 ft/sec 


= 176.468°. 


The equation set gives the patch point coordinates as 


6 nrt 


4.493° 


X NRT 


30.803°. 


Solving the same problem with a standard patched conic analysis, 
the patch point coordinates are found to be 


e NRT 


4.424° 


X NRT 


30.894° 


which agree well with the equation set calculated coordinates. 
Further demonstration of the validity of this equation set will 
be provided in Section 4.0. 
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As pointed out earlier, there is a lower limit on the 
energy of major-body-centered single-rectilinear trajectories 
and, consequently, on the energy for which this equation set 
will produce a solution. However, non-existence of these solu- 
tions does not mean that non-rectilinear trajectories do not 
exist at energies belov/ this level. Indeed, many trajectories 
of interest do exist with energy levels below the minimum 
energy of major-body-centered single-rectilinear trajectories. 

3 . 6 Non-Rectilinear Trajectory Patch Point Equations (General Set ) 

The energy limitation on the high energy equation set 
can be circumvented by eliminating the major-body-centered 
single-rectilinear trajectory patch point locus and deriving 
a new locus referenced to the minor-body-centered single- 
recti] inear trajectory patch point. Specifically, by assuming 
that the minor-body-centered trajectory plane passes through 
the minor- body-centered single-rectilinear trajectory patch 
point, it is possible to find an explicit statement for the 
patch point coordinates of a non-rectilinear trajectory. Since 
wc already have the equations of the minor-body-centered single- 
rectilinear patch point locus, we need only obtain one more 
equation set. Figure 7 shows the trajectory geometry required 
for non-rectilinear trajectories. From the geometry of the 
figure, the minor-body-centered state vector at the sphere of 
influence can be written as 


R, = R, 


COS SRT2 (COS6 SRT2 COSy NRT " Sin6 SRT2 Slni l SinY NRT* 


+ silU SRT2 COS1 l SinY NRT 


Sin> SRT2 (cOS6 SRT2 COSY NRT " Sin6 SRT2 Sini l SinY NRT ) 


- cosX SRT2 cosi lS inY NRT 


Sin£ SRT2 COSY NRT + COS6 SRT2 Sini l SinY NRT 


(53) 




t> 
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= V, 


I 

I c °SAgRT2 ( cos ‘ i SRT2 COS ’ NRT " sin3 SRT2 sini l SinY NRT* 

I 

+ sinA SpT 2 cosi 1 sin', NRT ] (-cos^) 

+ [cosx SRT2 (cos? SRT2 sin> NRT + sinL 'sRT2 sini l COSY N 

- sinA SRT 2 cosi^cosY NRT ] (sin<j> 1 ) 

[sinX SRT2 (COSf SRT2 COSY NRT _ Sin6 SRT2 sini l SinY NRT ) 

I 

- cosA SRT 2 cosi 1 sinY NRT ] (-cos*-,) 

+ [sinX SRT2 (c ° sB SRT2 SinY NRT + sin6 SRT2 sini l COSY N 

I 

+ cos.\ SRT 2 cosi 1 cosY NRT ] (sin^) 

lsin6 SRT2 COS 'NRT + COS D SRT2 s ^ n *l sinif NRT^ ( ' COS H ) 

+ [sin6 SRT2 Sin,, NRT " cos6 SRT2 slni i COS, NRT I (sin * 


where, from Figure 7, 


cosi. 


cosi-, = 


1 cos 6 


SRT2 


For this case, we must have 


R 2 x V 2 h 2 • 


The major-body-centered angular momentum vector may be written 



~sini 2 sinn NRT2 ' 

-sini 2 cosn NRT2 

_cosi 2 


(54) 


(55) 


(56) 


(57) 
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The minor-body-centered angular momentum vector, 
expressed in the major body coordinate system, is 


R]. x V]_ = h i 


-sini lS inO NRTi 


+sini 1 cosO NRT1 


cosi. 


(58) 


The angle n NRTl may be calculated from the minor-body- 
centered single-rectilinear patch point as follows : 


Combining Eqns. 
yields 


fi NRTl “ X SRT2 
(4) through (7) 



and Eqns. (53) through (58) 


(59) 




TSOO 



Bellcomm. Inc. 


32 


Equation set (60) provides three equations with the 
unknown quantities y nrt and 0 NRT2 - 

The third component of equation set (60) is a func- 
tion of > only. Moreover, it may be written in the form 

NKl 


Asin ’' N RT + BCOSY NRT + C * °‘ 


(61) 


Then i NRT is found, from Eqn. (13), to be 


siny 


-AC ± B 


NRT 




A 2 + B 2 - C 2 


”2 i 
A^ + B 


where 


(62) 


A = R 1 V o I( s in0 o )(cosX SRT2 sin6 SRT2 sini 1 - sinXg^cosi^ 

+ (cos* Q ) (sinX SRT2 sin$ SRT2 sini’ + cosX SRT2 cosi ’ ) ] 

+ Vl I(sin V (_sinX SRT2 COSB SRT2 ) 

I I 

+ (cos^q^) ( — sinX g RT2 s in 8 SRT2 sini ^ - cosXg^jCOsi^ 1 

B = R 1 V 0 l-sin4> 0 cosX SRT2 cos6 SRT2 - cos$ 0 sinA SRT2 cosB SRT2 ] 

+ RqV^ [ (sin^) (-sinXg^jSinBg^jSini^ - cosX sRT2 ;;os *l^ 
+ cost^sinX SRT2 cos 8 srt2 1 

C = h^cosi^ - h 2 cosi 2 + R 0 V 0 sin4> 0 


I 

cosi^ 


cosi^ 

cos6 SRT2 


and from Eqn. (55) 
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The quantity is found from 


hi = RiViSin^i* 


In these equations ii and ii are signed according to the rule: 

I 

il and ii are positive if the nearest node is descending, 

I 

il and ii are negative if the nearest node is ascending. 

Equation (62) provides two pairs of values of y nrt » the 

principal value and its supplement, for each choice of the sign on 
the radical. However, the supplementary values can be eliminated 
on geometrical grounds, since values of y nrt in excess of V0° 

would place the patch point in the wrong hemisphere. Selection 
of the correct value of y nrt from the two principal value solu- 
tions presents a more difficult problem. Figure 8 shows the 
varia'v ! n of y nrt with ii for a typical set of a^, h^, i 2 and h 2 

values xn the Earth-Moon system. Values for sinn NRT , and cosfi NRT2 

are obtained from the first and second components of equation set 
(60). Figure 9 shows the variation of the root sum square of the 
sine and cosine for the same parameter set used in Figure 8. It 
is seen that, for almost every value of i^, only one value of 

. NRT satisfies the complete equation set, that is, only one value 

°f > NRT yields a value of 1 for the root sum of the squares oi 

the sine and cosine of & NRT 2* 

However, Figure 8 shows two regions of i^ values for 
which no y nrt value is available, while Figure 9 shows six dis- 
crete points for which both values of y nrt satisfy the equation 
set. The second of these difficulties will be dealt with first. 


For a given set of a^ i 2 , and h 2 values, there are 

four possible minor-bodied-centered single-rectilinear patch 
points. These are shown in Figure 10. Tney are distinguished 
by the sign associated with i 2 and by the sign chosen on the 



1 m m m 


c 


FIGURE 8 - THE VARIATION OF y Hm WITH MINOR -BODY-CENTERED INCLINATION FOR A TYPICAL EARTH-MOON TRAJECTORY 
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FIGURE 10 - THE PATCH POINT GEOMETRY REQUIRED FOR A DOUBLE SOLUTION TO 7 NRT 
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radical in Eqn. (43) for x gRT 2 (see discussiv following 

Eqn. (45)). However , these distinctions are lost in the opera- 
tions used to solve equation set (60). Now, three minor-body- 
centered planes may be drawn joining the primary (i.e., intended) 
and secondary minor-body-centered single-rectilinear patch 
points, also shown in Figure 10. When i^ matches the inclina- 
tion of these planes, two values of y nrt are possible and these 

are the two values found by Eqn. (60). There are six points with 
double solutions since each of the three planes may be traveled 
in either a retrograde or posigrade sense. 

The other anomaly in Eqn. (62) is the absence of any 
value of y NRT for certain values of i^. Specifically, the 

radical in Eqn. (62) becomes imaginary for these values of i^. 

This occurs because the equations are based on the assumption 
that the minor-body-centered plane passes through the minor- 
body-centered single-rectilinear patch point. While this is a 
good approximation, it is not strictly true. As a result, in 
the region where the radical should be nearly zero, it becomes 
imaginary for a small range of i^ values. 

We now have, with the exceptions noted, a unique 
solution for y^ r <t* Using basic spherical trigonometry 

relationships , 


sinA6 NRT = sini l sin YNRT (63) 

and 

I 

tanAX NRT = cosi i tan YuRT * (64) 

The non-rectilinear patch point coordinates are then 


and 


®NRT * e SRT2 + AS NRT 

► 

X NRT ” X SRT2 " AX NRT* 


(65) 
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As in the high energy equation set, an alternative 
formulation for the non-rectilinear patch point may be obtained 
by calculating the angular distance from the minor-body-centered 
node line to the minor-body-centered, single-rectilinear patch 
point, adding to this y nrt and then calculating the non- 


rectilinear patch point coordinates using this sum and i^. The 
resulting equations are 


Sin8 NRT = sini l sin 

''NRT + Sin ’ 1 1 

fsingsRT 2 \ 
i sinl i | 





and 


> ( 66 ) 


NRT 


= X SRT2 + Sin 


-1 


/tang 


SRT2 


tanx. 


sin 


-1 




j 


As there seems to be no advantage in accuracy to one formulation 
over the other, the simpler Eqns. (65) are used here. 


Figure 11 is a plot of Eqns. (65) , with the locus for 
a formal patched conic trajectory included for comparison. The 
loci are in excellent agreement except in the region where y NRT 

does not exist and its immediate neighborhood. In this neighbor- 
hood the locus calculated by Eqns. (65) deviates from the for- 
mally calculated patched conic locus. 


From a practical calculative viewpoint, it is desirable 
to prevent the loci from diverging and to obtain a value for 
, NRT in the imaginary region. An algorithm is derived to accom- 
plish this purpose, based on the following observations. First, 
the real locus is nearly symmetrical about a great circle line 


of inclination 8 SRT2 


and passing through the minor-body-centered 


single-rectilinear trajectory patch point. In fact, for all but 
very low energies, it is nearly circular.* At the same time, 
when the function generated locus is misbehaving for a given 
value of i 1# it is usually well behaved for -i^. Thus, using 


the symmetry of the locus, one may obtain a solution for Y NRT 


by changing the sign on i^. 


* Indeed, the major-body-centered single-rectilinear trajec- 
tory locus (Eqn. (23)), used in the high energy equation set to 
approximate the non-rectilinear locus, is circular, since it is 
centered at the dual rectilinear patch point and is not a function 
of i r 
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The next problem is deciding when the function is 

2 2 2 

misbehaving. This always occurs when A + B - C is small 
2 

compared to C . Thus, one may make the rule that i, is set 

2 2 2 2 

equal to -i^ in Eqn. (73) if A + B - C < eC where e is 

an empirically derived small number. The author has used 
e = .01 with success. However, this rule is not quite adequate, 

2 2 2 2 

as the function may yield values of A + B - C < eC for both 
+i 1 and -i^; also, it may generate imaginary values for 




+ B - C for both i.'s. These problems can be handled as 


follows: if the radical is imaginary for both +i^ 'nd -i^, 

set it equal to zero; if only one value of the rar ,il is 
imaginary, use the value of i^ that yields the rec answ' r; 

and if both values of i, yield real answers but still 

<"» a 


2 2 
h + B 


2 2 

- C < eC , use the i, that produces the larger 
A 2 t B 2 - C 2 . 


Stated in algorithmic form, the above rules are: 
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Note that the sign of ^ is changed in Eqn. (62) only. Equa- 
tions (63) and (64) use the proper value of i^. 

It is now possible to select the proper value of y nrt 

from the two solutions of Eqn. (62) without resorting to the 
JL NRT2 ec 5 uat i° n s. Observe from Figures 8 and 9 that the proper 

value of y nrt is small compared with the extraneous solution. 

The only region where this is violated is the neighborhood of 
the imaginary region, but this region has been eliminated by 
the above algorithm. Consequently, we may simply state: the 

sign on the radical in Eqn. (62) is chosen to obtain the smallest 
positive value of y nrt for trajectories from the minor body to 

the major body, and the smallest negative value for trajectories 
from the major body to the minor body. 

Figure 12 presents the same locus shown in Figure 11, 
incorporating the corrective algorithm given above. Also shown 
is the true patched conic locus. As can be seen, the locus is 
prevented from deviating seriously away from the path of the 
real locus. 


The difficulties just discussed in the solution of 
Eqns. (70) become worse as the energy is decreased. To demon- 
strate that the equation set plus algorithm operates satis- 
factorily at very low energies, Figure 13 shows the equation 
set locus and the true patched conic locus for a low energy 
(116 hour) trip to the Moon. Comparison with the true locus 
shows the equation set to be producing quite good answers. 


Equation (62) is limited in energy range by the minimum 
value of the minor-body-centered single-rectilinear trajectory 
energy. That is, y nrt cannot be calculated from Eqn. (62) for 


energies below the minimum for which a value of 


X SRT2 exists ' 


as given in Eqn. (45) . Do such low energy trajectories exist? 

If so, how large is the band of energy missed by Eqn. (62)? 

To answer the first question, consider the nature of the deriva- 
tion. Only one assumption was made, i.e., the non-rec til inear 
trajectory plane centered at the minor body passes through the 
minor-body-centered single- rectilinear patch point (& SRT 2' X SRT2^ 


If this assumption were strictly true, then the limiting energy 
for * srt2 would also be the limiting energy for non-rectilinear 


trajectories. However, this is not the case. 




aaruuvr 


FIGURE 12 - PATCH POINT LOCUS GENERATED BY THE GENERAL EQUATION SET 
USING THE CORRECTIVE ALGORITHM 




CONIC 
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Observe the case of major-body-centered, single-rectilinear 
trajectories. Here the associated value of Ag RT2 must 

dual rectilinear patch point. For this case, it was shown that 
the minor-body- centered plane passed through the dual recti- 
linear patch point, but the proof required certain approxima- 
tions; thus, the result is only approximately true. Conse- 
quently, non-rectilinear trajectories may exist with energies 
below the minimum major-body-centered single-rectilinear 
trajectory energy; in fact, such trajectories do exist. 


How wide a band of energy is unavailable to Eqn. (62)? 
Since the general equation set yields accurate results for quite 
low energies as demonstrated by Figure 13 , one may argue that 
the assumption used must be fairly accurate and thus the missed 
energy band is quite narrow. Consider the minimum energy tra- 
jectory to the Moon (at apogee) which the general set will cal- 
culate. Take as input: 


i 2 » 30° 

h 2 = 7.72 x 10 11 f t 2 /sec 
i x = 160° 

R _ = 6.06698 x 10 6 ft (60 n.mi. altitude). 

r D 

Then the minimum energy direct trajectory has 
a x = 34.8609 x 10 6 ft 

and a trip time of 134 hours. Based on Apollo trajectory studies; 
direct trip times noticeably in excess of 134 hours do not exist. 
Thus, it is safe to say that the band of trajectory energies 
missed by Eqn. (62) is very narrow. 

As in the high energy equation set, there are six 
trajectory elements, or independent parameters, which deter- 
mine the spacecraft trajectory. They are: 

minor-body-centered 

inclination (i^) 

velocity at the sphere of influence (V^) 
trajectory flight path 


angle t^) 


\ 
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ma j or-body-centered 

inclination (ij) 
angular momentum (h 2 ) 

and the time of piercing the sphere of influence. 


The elements V^ and may be restated in terms of 

other parameters, if desired, using the standard conic rela- 
tionships. For example, 



(67) 


^1 


T1 


♦ 2 


sin 


^1 

— + 2 
a i 


( 68 ) 


obtains and from the semi-major axis and radius of peri- 

center. Means of obtaining h 2 are discussed in Appendix C. 

R q , V , and $ 0 are determined from the minor-body ephemeris 

and calculated at the time the spacecraft pierces the sphere 
of influence. 

The general equation set is valid over essentially 
the entire regime of possible energies. While there is a limit 
on the validity of the . equation, it is seen that this is 
not significant. 

4.0 VALIDITY OF THE PATCH POINT EQUATIONS 

In order to demonstrate that the techniques presented 
here are a good approximation to a rigorous patched conic 
analysis, a series of trajectory problems have been solved using 
a patched conic analysis, the high energy equation set, and the 
general equation set. The results are shown in Tables I through 
VI. Each table presents the results for three or four different 
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COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND THE HIGH ENERGY 
AND GENERAL EQUATION SETS FOR HIGH ENERGY LUNAR TRAJECTORIES 



SPHERE OF INFLUENCE 
PA'i CH POINT 


PERI- 

LATITUDE LONGITUDE CENTFR 
RADIUS 


MINOR BODY 
PARAMETERS 


INCLI- 

NATION 


MAJOR BODY 
PARAMETERS 


SEMI- 

MAJOR 

AXIS 


DEGREES DEGREES FEETxlO 6 DEGREES FEET x 10 


PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 



26.000 6.0670 

25.992 6.0670 

25.006 6.0670 



PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 



160.012 

160.012 

160.012 


PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 



... 1 — 

26.145 

26.145 

26.145 

6, 

6, 

6, 




PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 



134.841 

134.841 

134.841 


EARTH-MOON TRA JECTORIES WITH A TRIP TIME OF APPROXIMATELY 60 HOURS. THE MOON IS AT 
APOGEE. THE TRAJECTORY PERIGEE ALTI i UDE IS 100 N.M. AND PERISELENE ALTITUDE IS 60 N.M. 



TABLE II 


COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND THE HIGH ENERGY 
AND GENERAL EQUATION SETS FOR INTERMEDIATE ENERGY LUNAR TRAJECTORIES 
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TABLE III 

1 COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND 

l THE HIGH ENERGY AND GENERAL EQUATION SETS FOR LOW ENERGY LUNAR TRAJECTORIES 



SPHERE OF INFLUENCE 
PATCH POINT 

MINOR BODY 
PARAMETERS 

MAJOR BODY 
PARAMETERS 

LATITUDE 

LONGITUDE 

PERI- 

CENTER 

RADIUS 

INCLI- 

NATION 

SEMI- 

MAJOR 

AXIS 


INCLI- 

NATION 

DEGREES 

DEGREES 

FEET x 10 6 

DEGREES 

FEET x |0 6 

FT 2 /SECx10 U 

DEGREES 

PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 


58.500 

6.0670 

44.102 

40 SOLUTlOf 
44.102 

22.935 

i 

1551 

23.119 

8.173 

58.571 

6.0670 

i - 

22.935 

1119 

23.879 









PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 

8.000 

74.500 


159.324 
JO SOLUTlOf 
159.324 

23.718 

■ 

7.7175 

31.980 

8.475 

74.766 


i 

23.718 

7.8933 

33.857 








■■■■■■■ 

PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 

8.000 

7.879 

7.845 

66.000 

65.880 

65.954 

6.0670 

6.0670 

6.0670 

25.367 

25.367 

25.367 

in 

7.7175 

7.6357 

7.6277 

135.072 

135.554 

135.803 








■■■■ 

PATCH CONIC 

HIGH ENERGY EQN. SET 

GENERAL EQN. SET 

8.000 

7.931 

7.882 

77.000 

77.204 

77.135 

6.0670 

6.0670 

6.0670 

143.992 

143.992 

143.992 

11.608 

11.608 

11.608 

7.7175 

7.7040 

7.6742 

140.265 

140.760 

140.982 


EARTH-MOON TRAJECTORIES WITH A TRIP TIME OF APPROXIMATELY 90 HOURS. THE MOON IS AT PERIGEE. 
THE TRAJECTORY PERIGEE ALTITUDE IS 100 N.M. AND PERISELENE ALTITUDE IS 60 N.M. 


TABLE IV 

COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND THE 
HIGH ENERGY AND GENERAL EQUATION SETS FOR MARS TRAJECTORIES 





























































































-48- 


I 

TABLE V 


COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND 
THE HIGH ENERGY AND GENERAL EQUATION SETS FOR JUPITER TRAJECTORIES 



SPHERE OF INFLUENCE 
PATCH POINT 

MINOR BODY 
PARAMETERS 

MAJOR BODY 
PARAMETERS 

LATITUDE 

LONGITUDE 

PERI- 

CENTER 

RADIUS 

INCLI- 

NATION 

SEMI- 

MAJOR 

AXIS 


INCLI- 

NATION 

DEGREES 

DEGREES 

A. U. 

DEGREES 

A. U. 

(A.U.) 2 /DAY 

DEGREES 

PATCH CONIC 

HIGH ENERGY EON. SET 

GENERAL EON. SET 

0.800 

-64.000 

.00465 

1 

33.403 

MO SQLUTIOr 
| 33.403 | 

0.033897 

4 

0.022370 

0.559 

1.233 

-63.861 

.00465 j 

| 0.033897 


0.234 



■■BBBI^HB 



BI^^BBBBI 

^BHBBBBB 

PATCH CONIC 

HIGH ENERGY EON. SET 

GENERAL EON. SET 

0.800 

-74.000 

.00465 

f 

| 134.183 | 0.035344 

an cm i iTiniu 

0.022370 

0.935 

1.628 

-73.725 

.00465 | 

| 134.183 | 

0.035344 

0.022367 

0.316 








PATCH CONIC 

HIGH ENERGY EON. SET 

GENERAL EQN.SET 

1.600 

-66.000 

.00465 | 

_ * 

| 60.829 | 

JO SOLUTIOI^ 
| 60.829 | 

| 0.034210 

j 

0.022370 

0.641 

2.347 

-65.843 

.00465 j 

0.034210 

0.022369 

0.159 




mi 





PATCH CONIC 

HIGH ENERGY EON. SET 

GENERAL EQN. SET 

0.600 

0.634 

0.846 

-86.000 

-86.402 

-85.887 

.00465 

.00465 

.00465 

138.163 

138.163 

138.163 

0.0022498 

0.0022498 

0.0022498 

0.022370 

0.022390 

0.022388 

176.334 

176.242 

175.659 


EARTH-JUPITER TRAJECTORIES WITH A TRIP TIME OF APPROXIMATELY 945 DAYS. JUPITER IS AT APHELION, THE 
TRAJECTORY PERIHELION RADIUS IS 1 A.U. AND THE PERICENTER RADIUS AT JUPI TER IS 500,000 N M. 


f 

i 

I 

TABLE VI 


COMPARISON OF RESULTS FROM PATCH CONIC ANALYSIS AND THE 
HIGH ENERGY AND GENERAL EQUATION SETS FOR JUPITER TRAJECTORIES 



SPHERE OF INFLUENCE 
P*TCH POINT 

MINOR BODY 
PARAMETERS 

MmJOR BODY 
PARAMETERS 

LATITUDE 

LONGITUDE 

PERI- 

CENTER 

RADIUS 

INCLI- 

NATION 

SEMI- 

MAJOR 

AXIS 


INCLI- 

NATION 

DEGREES 

DEGREES 

A.U. 

DEGREES 

A.U. 

(A.U.) 2 /DAY 

DEGREES 

PATCH CONIC 

HIGH ENERGY EON. SET 

GENERAL EQN. SET 


-66.000 

.00465 

26.141 

JO SOLUTIOI 
26.141 

0.034783 

.022366 

0.333 

1.214 

-65.838 

.00466 

0.034783 

.022366 

0.078 


■bbbbb 

^hbi 

h^bi^hbbbi 

BBBBBBBBi 

BBBBB^^BI 

BBRBBBBBI 

bbbbbbb. 

PATCH CONIC 

HIGH ENERGY EQN.SET 

GENERAL EQN.SET 

1 

-66.000 

.00465 

25.949 

JO SOLUTIOI 
25.949 

0.034727 

.022360 

1.242 

0.040 

-65.800 

.00466 

0.034929 

.022364 

.906 







■■■■■■■■■■■ 

^BBBBBBBM 

PATCH CONIC 

HIGH ENERGY EQN.SET 

GENERAL EQN. SET 

0.800 

-82.000 

.00466 

f 

131.892 
K) SOLUTION 
131.892 

0.036691 

i 

.022366 

1.015 

3.363 

-80.916 

.00465 

0.036691 
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trajectories, and for each trajectory, three rows of values 
are presented. The first row contains the results of the 
patched- conic analysis. These results constitute the input 
values for the high energy and general equation sets. The 
patch points calculated by the equation sets are given in the 
first and second columns of the last two rows in each table. 

These patch points should (and do) agree closely with the 
patched- conic patch point in the first row. Finally, the state 
vector obtained using the input minor-body-centered parameters 
and the patch points calculated by the equation sets were used 
to calculate the major-body-centered inclination and angular 
momentum. These values are presented in the last two columns 
of the second and third rows of each table. Again, these values 
should (and do) agree well with the patched conic angular mo- 
mentum and inclination of the first row. The trajectory prob- 
lems solved cover a range of Earth-Moon, Earth-Mars, and Earth- 
Jupiter trajectories. For the interplanetary trajectories, the 
Earth-centered trajectory portion was ignored and the trajectories 
were assumed to have a 1 A.U. perihelion radius. 

Tables I through III cover Earth-Moon trajectories for 
a range of geocentric energies. The sphere of influence radius 

O 

was taken as 1.8 x 10 feet for consistency with BCMASP.* Study 
of these tables shows that the high energy equation set produces 
excellent results in those cases where it is valid. However, 
it is not capable of calculating the patch point for low energy 
trajectories within the range of interest. 

The results generated by the general equation set are 
comparable in accuracy to those from the high energy equation 
set. Moreover, the general set is able to calculate nearly all 
the trajectoi ies unavailable from the high energy set. 

Mars and Jupiter trajectory results are tabulated in 
Tables IV, V, and VI, For this set the radius of the sphere of 
influence was calculated as 



( M minor body ! R 
M major body] 1 


(69) 


The results for Mars (Table IV) were very good. The excellence 
of these results comes from Mars' small sphere of influence and 
the small pericentar distance chosen for the problem. The re- 
sults for Jupiter (Table V) are not quite as accurate. Jupiter 


*Bellcomm Apollo Simulation Program 
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has a large sphere of influence compared with Mars and the 
pericenter radius chosen for the study was also relatively 
large. Even so, the results are satisfactory. 

Table VI shows four more trajectories from Earth to 
Jupiter. These have a slightly longer trip time than those in 
Table V and are very nearly equivalent to a Hohmann transfer. 
The results of the first three trajectories, which have posi- 
grade helicentric portions, are similar to those of Table V. 
The last trajectory, which has a retrograde heliocentric por- 
tion, could not be calculated using either equation set. 


The failure of both equation sets to calculate a patch 
point for this trajectory is due to the limits on the minor- 
body-centered single-rectilinear trajectory energy noted in 
Eqn. (45). The minor-body-centered energy selected is too 
small to give a real solution for ^ SRT ^ i n Eqn. (43) . 


To find the real lower limit on the Jupiter trajectory 
energy imposed by Eqn. (45) , the input value of a^ was varied 

to find its minimum value for an Earth- Jupiter trajectory. The 
other trajectory elements selected were: 


Rpi = 500,000 n. mi. 
i 2 = 1.3“ 

h 2 = .02236 (A.U. ) 2 /day. 


These i lements might represent a mission to Jupiter's moon 
Europa (Jupiter II) . Jupiter was taken to be at aphelion. The 
minimum value for a^ was found to be .03688 A.U. This repre- 
sents a trajectory aphelion distance of 5.438 A.U. compared to 
Jupiter's heliocentric distance of 5.455 A.U. Thus, the helio- 
centric energy found is actually slightly less than the energy 
of a formal Hohmann transfer, so we can say that the general 
equation set does indeed cover essentially the entire spectrum 
of possible Earth- Jupiter trajectory energies. 

5.0 SUMMARY 

Two equation sets have been derived which provide a 
good approximation to the problem of patched conic analysis of 
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spacecraft trajectories between two large central bodies satis* 
fying boundary conditions imposed at both central bodies. A 
general equation set provides solutions over essentially the 
entire spectrum of interesting trajectories, failing only for 
extremely low energy trajectories. A high energy equation set 
provides solutions only for the higher end of the trajectory 
energy spectrum; however, it is somewhat simpler than the 
general set. For easy reference, these equation sets are or- 
ganized into a computational algorithm in Appendix B. 

The primary benefit obtained from these equation sets 
is an exact statement of the influence of the different vari- 
ables on the trajectory problem. The analyst wishing to come 
to grips with any problem of trajectory design is provided with 
a powerful tool. By using derivatives of the patch point loca- 
tion, problems in post-periselene abort and midcourse correc- 
tion should become much more tractable. The problem of free 
return and planetary swingby trajectories can also be reduced 
to analytical form (Reference 2). 

Another significant benefit is the improved calcula- 
tion time for patched conic trajectories. Not only is the time 
on a computer reduced, but the non-iterative solution format is 
appropriate for use on a desk calculator if desired. 

9f $*J*4*y 

2011-KMC-vh K. M. Carlson 

Attachments 

Appendixes A, B, and C. 
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APPENDIX A 


DEFINITIONS , SYMBOLS, AND SUBSCRIPTS 

A number of terms, symbols and subscripts are 
utilized in the text. Definitions of frequently used terms, 
as well as a listing of symbols and subscripts used, are 
repeated here for easy reference. 

A. 1 Definitions 

Major Body - The more massive of the two large central 

bodies. The gravitational field of this body 
is assumed to pervade all space except a spheri- 
cal region about the minor body. 

Minor Body - The less massive of the two large central bodies. 

The gravitational field of this body is assumed 
to act only inside its sphere of influence , a 
spherical region of space centered at the minor 
body. The radius of the sphere of influence is 
taken to be 

„ _ minor body! 2/ ^ 5 „ 

R 1 Iff . . I R o’ 

\ major body/ 

A trajectory which rises and/or descends verti- 
cally with respect to its central body. Such 
trajectories have zero angular momentum. 

A trajectory between two central bodies in 
which the portions of the trajectory centered 
at each body are rectilinear. 

A trajectory between two central bodies with 
only one rectilinear body centered portion. 


A trajectory which has no rectilinear portion. 
This is, of course, the most general case. 


Rectilinear - 
Trajectory 


Dual 

Rectilinear 

Trajectory 

Single 

Rectilinear 

Trajectory 

Non- 

Ractilinear 

Trajectory 
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A. 2 Symbols 

a “v semi-major axis 
h angular momentum 
i ^ inclination 

, , cosi, 

i, ^ modified inclination, i, = —— e , taken to have 

1 1 COS0 SRT2 
the same sign as i^ 

J+l for trajectories from the minor body to the major body 
k W 

1-1 for trajectories from the major body to the minor body 
M ^ mass 
R ^ radius 
V ^ velocity 

2 ^ latitude 

y ^ the distance the patch point is displaced from its 
reference point 

A 'v. longitude 

v ^ gravitational constant of the body 

♦ ^ the angle between the radius vector and the velocity 
vector , measured from the radius vector 

ft ^ the angle between the trajectory node nearest the 

sphere of influence pierce point and the line between 
the major and minor body centers 

A. 3 Subscripts 

DRT % dual rectilinear trajectory 
SRT1 'v major-body-centered single-rectilinear trajectory 
SRT2 -v minor-body-centered single- recti linear trajectory 
NRT -v non-rectilinear trajectory 
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0 'v a minor body property referenced to the major body 

1 ^ a spacecraft property referenced to the minor body 

and occurring at the sphere of influence 

2 'v a spacecraft property referenced to the major body 

and occurring at the sphere ' f influence 

E ^ Earth centered 
S ■v selenocentric 
A ^ apocenter 
P ^ psricenter 
AG ^ apogee 
PG - perigee 
PS •• periselene 
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A COMPUTATIONAL ALGORITHM FOR PATCHED CONIC TRAJECTORIES 

The equation sets have Leen organized into an algo 
rithmic form to facilitate their use. The required inputs 
are listed and the calculation described in a step by step 
fashion until the final solution is determined. 

B. 1 Basic Data Required 

u 1 » y 2 " the gravitational constants of the major and 
minor bodies 

R q - the distance between the major and minor body 
at the time of penetration of the sphere of 
influence 

$ - the angle from the minor body radius vector to 

the minor body velocity vector at the time of 
penetration of the sphere of influence 

V Q - the velocity of the minor body at the time of 
penetration of the sphere of influence 

B.2 Trajectory Inputs 

- the spacecraft minor-body-centered velocity at 
the sphere of influence 

- the spacecraft minor-bcdy-centered flight path 
angle of the sphere of influence 

i. - the spacecraft minor-body-centered inclination, 
referenced to th,i minor body's orbital plane 

h- - the spacecraft major-body-centered angular 
momentum 

i 2 - the spacecraft major-body-centered inclination, 
referenced to the minor body's orbital plane 

B. 3 Alternative Inputs 

- the minor-body-centered semi-major axis? V, is 
obtained from 
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V 


1 


2 



R pl - the minor-body-centered radius of pericentron; ^ 
is obtained from 


. 2 
sin 





R A2 and Rp2 - the major-body-ce^ered apsidal distances; 

is obtained from 


h 


2 



u 2 R A2 R P2 
R A2 + R P2 


Note that if R_ 0 << R.,, as is typical of Earth-Moon trajectories, 
then 


h 2 ^ *^ W 2 R P 2 


The sphere of influence radius is obtained from 
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B * 4 The Minor-Body-Centered Single-Rectilinear Patch Point 


Sine SRT2 = - 


h 2 sini 2 sinij SRT ., 

R l V o sln *o 


where 


tanft 


SRT2 


R l V o cos *o + R o kV l 


and 


smx 


-AC ± B 


SRT2 


bV 


A 2 + 


B - C" 




B' 


where 


A = R kV, + R..V cosd) 

0 1 1 o Y o 

B = R.V sin<f> 

1 o y o 


C = - 


R o V o sin *o " h 2 cosi 2 


cose 


SRT2 


The angles 6 g RT 2 and *g RT 2 are roinor-body-centered single- 

rectilinear patch point coordinates. The sign in the x-__- 
, , SRT2 

equation is chosen according to 


a) 

b) 


X SRT2 < sin 1 (”f| ' 
X SRT2 > sin 1 ("?) ' 


use the minus sign, 
use the plus sign . 
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The choice of sign differentiates the minor-body-centered 
single- rectilinear trajectories roughly into those lying com- 
pletely within the minor body's orbital path (Case a) and 
those that lie partly or completely outside the minor body's 


trajectory path (Case b) . 


SRT2 


lies in the first or second 


quadrant for trajectories from the minor body to the major 
body, and in the third or fourth quadrant for trajectories 
from the major body to the minor body. 


B. 5 The Patch Point Displacement 

Next, calculate the trajectory patch point displace- 
ment from the minor-body-centered single-rectilinear patch point 
by either of the following two methods. 


B.5.1 High Energy Trajectory Case 


siny 


V^sin$^ 


SRT1 


V sin* SinX DRT 
o Y o 


and 


sim D R T 


-AC ± B 


bV a 2 + 


2 2 
B Z - C^ 


1 T 

A + B 


f* 


with 


A = R 0 kV 1 + R 1 V 0 COS* 0 

B = R 1 V 0 sin^ 0 

C * -R Vsin<j> 
o o T o 

where k determines the direction of and obeys the rules: 

k = +1 for trajectories from the minor body to the major 
body 

k = -1 for trajectories to the minor body from the major 
body. 
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The angle ^ drt t * ie l° n 9 itude of the associated dual recti 

linear trajectory. The sign in the a drt equation is chosen 
according to UK1 


a) 

X DRT < S ^ n 

l-ll 

, use the minus 

sign 

b) 

if x DRT > sin -1 

l-ll 

, use the plus 

sign 


This procedure will produce values for y srt ^ only 

for relatively high energy trajectories. The following alter- 
native procedure is more general. 

B.5.2 General Trajectory Case 


siny 


NRT 


-AC 


. bV 


A 2 ♦ 


A + B 





with 


A = 


R l V of ai “*0 ) * COsX SRT2 sin8 SRT2 sini l ” sinX SRT2 COsi l^ 

+ (cos* G ) (sinA SRT2 sine SRT2 sini^ + cosA SRT2 cosi|)j 

+ R q V 1 j(sin<t> x ) (-sinA SRT2 cosB SRT2 ) 

+ (cos^) (-sinA SRT2 sinB SRT2 sini| - cosX SRT2 cosi^)J 

B - R 1 V C [-sin* 0 cosx SRT2 co S B SRT2 - cos* 0 sinx SRT2 co.B SRT2 ] 

+ R 0 V 1 [(sin^l (-sinx SRT2 sine SRT2 «inij - cosx SRT2 co»i’ ) 

_l_ At.... . — J u . AAA . "l 


c = 


+ cos* 1 sinx SRT2 cosB SRT2 

m 

I 

h^cosi^ - h 2 cosi 2 + R 0 V 0 sin$ 0 


where 


cosi. 


cosi. 


CO>S SRT2 ' 
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and 


= R 1 V 1 sin<(>^ . 

I 

The quantities i^ , i^ and i 2 carry signs, according to the rules 


— 

Nearest minor-body-centered node 

Ascending 

Descending 

• 

i i and i^ are negative 

1 

i^ and i^ are positive 


and 


Nearest major-body-centered node 

Ascending 

Descending 

i 2 is positive 

i 2 is negative 


Having chosen the appropriate sign for and 

by the above rule, this sign must be modified (for both i, and 
• 1 

i^ ) for use in the y nrt equation according to the following 
algorithm: 
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Kote that the sign of is changed in the y nrt 

equation only. The sign on the radical in the y nrt equation 

is chosen to obtain the smallest positive value of y nrt for 

trajectories from the minor body to the major body and to 
obtain the smallest negative value of y nrt for trajectories 

from the major body to the minor body. 

B. 6 The Non-Rectilinear Patch Point 

The non-rectilinear patch point may be obtained from 
the values obtained so far in either of two ways: 
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a) The patch point latitude is given by 

0 NRT = 6 SRT2 + a6 

and the longitude by 


A NRT A SRT2 " AA * 


where 


sinas * sini i sin v SRT1 


tan^A = cosij^ tany SRTl 


or 


sinA6 » sini x siny NRT 


tanAA = cosi^ tany NRT 


depending on whether Yg R »p]_ or y^ R ,p was obtained in the previous 
stsp* As bsforG 


COSl, 


COS1-1 = 


1 cos 6 


SRT2 


b) The alternative formulation for A NRT is 


s ^ n ^NRT s Sinitsin 


V SRT1 + sin 


-1 f sine SRT 2l 
sini. 
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or 


sine NRT = sini^sm 


y N rt + sin 


-1 ( Sin6 SRT2\ 


sini. 


depending on whether y SRT i or y nrt was obtained in the previous 
step. Finally 


. _i_-l | tan6 SRT2l 
X NRT X SRT2 + Sin tan i. 


- sin 


•1 l tan8 HRT\ 
\ tan 1 l) 


B. 7 The State Vectors at Sphere of Influence Penetration 
B.7.1 The Minor-Body-Centered State Vector 

The longitude of the nearest node line is 


. .< 'I ‘““hRTI 
n l " X NRT + • ln tan i, 


The argument of latitude of the patch point with respect to 
is 


-i sinS NRT\ 

“1 = sin TTnT; 


The position vector is then 


X 1 = R i (cosn^cosai^ + sinh^cosi^sinai^) 
^1 “ R i (sinn^cosa)^ - cosh^cosi^sinu^) 


- R^sini^sinu^ 
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The velocity vector is 

= V^[sin<t>i (cosfi^sinw^ - sinn^cosi^cosa)^) 

- cos<t^ (cosn^cosa)^ + sinfJ^cosi^sinw^) ] 

• 

= V^[sin4>^ (sinn^sinw^ + cos^cosi^cosw^) 

-cos4>^ (sinn^cosw^ - cosft^cosi^sinto^) 1 

Zj_ = [-sin<t>^sini^cosu>^ - cos^^sini^sinu)^] 

B.7.2 The Major-Body-Centered State Vector 


- V 0 cos * 0 
+ V 0 sin * 0 



The solution is now complete 
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APPENDIX C 


SOME USEFUL TECHNIQUES AND APPROXIMATIONS 

C. 1 Use of the Equation Sets as a Rigorous Patched-Conic 
Trajectory Generator 

For studies that require more precise results than 
those available from the equation sets, these sets may be used 
in an iterative scheme to provide a rigorous patched- conic- 
result. This has the advantage over standard patched- coni, o 
analyses that the iterative parameters are true trajectory 
parameters rather than the patch point coordinates. 

To create a rigorous patched-conic result, the major- 
body-centered angular momentum and inclination used in the 
equation set become shaping parameters, and are varied until 
the angular momentum and inclination that are actually obtained 
match the desired values. A simple mirror imagining technique 
is effective. That is, the amount the resultant ij and hj 

values exceed the desired values is subtracted from the input ij 

and values. The iteration may be performed on ij and hj 

simultaneously. Convergence is very rapid. 

C. 2 Approximating the Major-Body-Centered Angular Momentum 

In general, trajectory angular momentum is not one of 
the basic design parameters for a mission. Therefore, it would 
be useful to find a relationship between angular momentum and 
some other parameter which is basic to the mission design. For 
many problems of interest, this is easily done. 

C.2.1 Earth-Moon Trajectories 

The geocentric angular momentum may be written 

- (C-l) 

a E I 


a E . This can be seen by 
and that R AG ^ R Q - R^. 


VUfR 


E PG 


Rp G is typically much smaller than 
considering that a E * y(Rp G + R AG ) 
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Thus, the minimum a E is j(R pG + R 0 ” R^). Typical 
values for these parameters are 

R pG = .215 * 10® ft (100 n. mi. altitude) 

g 

R = 11.7 * 10 ft (closest approach) 

R x = 2 x io 8 ft 

In this case R pG is only 4.3% of a E . We may then make the 
approximation 


h 2 " V 2u E R PG* 


(C-2) 


This is the expression for the angular momentum of a parabolic 
trajectory with perigee radius R pG « Thus, typical Earth-Moon 

trajectories are nearly parabolic. Normally, the perigee radius 
can be assumed equal to the parking orbit radius for trajecto- 
ries from the Earth to the Moon. For trajectories from the Moon 
to the Earth which must enter the Earth's atmosphere in the 
Apollo entry corridor, a perigee altitude of 100,000 feet has 
been found to produce the desired results. 

The validity of this approximation has been tested 
for a number of trajectories. The results are tabulated in 
Tables C-I to C-VIII and summarized in Table C-IX. 

The data in Tables C-I to C-VIII were generated in 
the following steps. The title in parentheses for each step 
corresponds to the row name in the tables. 

a) (Pacched Conic) Calculate a set of trajectories 
using a rigorous patched-conic analysis for given 
values of Rp S , Rp G , geocentric energy and sphere 

of influence patch point coordinates. 

b) (Equation Set) Calculate the patch point with the 
equations set, using the values of 
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a l ' i l' R PS' 1 2' R PG 
obtained in (a) and calculating h E by 


= V 2p E R PG 


c) (inverse) Use the values of geocentric energy and 

patch point calculated in (b) , along with the values 
of Rpg and R p( , used in (a) , to make another rigorous 

patched-conic analysis. 

a) (Error) Take the difference between the values 

obtained in (a) and those obtained in (b) and (c). 

Finally, the average error between (a) and (b) is calculated and 
stated at the bottom of each Table. The Moon's parameters were 
calculated from the assumed lunar elements: 

g 

semi-major axis = 12.54 * 10 feet 
eccentricity = .06380 

Step (c) is intended to give the analyst a feeling 
for the significance of the patch point errors. 

Table C-IX is a listing of the average and maximum 
errors for each of Tables C-I to C-VIII. The results show 
very good agreement between techniques. 

C.2.2 .interplanetary Trajectories 

Another expression for angular momentum is 

h - <c-3> 


where R A is the apocenter radius and R p is the pericenter radius. 
Thus, for interplanetary trajectories angular momentum is related 
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to the trajectory perihelion and aphelion distances. It is 
usually more meaningful to choose the apsidal distances rather 
than the angular momentum of the trajectory in mission analyses . 
Also, many interplanetary trajectories closely approximate the 
minimum energy Hohmann transfer, and for such trajectories, 
the heliocentric distances of the planet of departure and the 
planet of arrival become the apsidal distances. 



TABLE C l 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION 

HIGH ENERGY EQUATION SET 
MOON AT APOGEE; 80 HOUR TRIP TIME 


TO ANGULAR MOMENTUM 


Equation Set 

Inverse 

Error 


Patched Conic 
Equation Set 


Patched Conic 
Equation Set 



Patched Conic 
Equation Set 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 


Patched Conk: 
Equation Set 


Patched Conk 
Equation Sat 


Patched Conk 
Equation Set 


0.00128 x 10 


0.14587 x 10 


0.14388 x 10 
0.00211 x 10 1 


0.15376 x 10' 


0.15171 x 10° 
0.00205x10* 


0.15571 x 10* 

0.15333 x 10* 
0.00235x10* 


GEOCENTRIC 

ENERGY 


-8.1700 x 10 6 
-8.1625 x 10 6 


-8.1700 x 10° 
-8.1609 X 10 6 


0.0091 x 10° 


-8.1700 x 10° 
-8.1591 x 10* 


0.0109 x 10° 


-8.1700 x 10° 
-8.1606 x 10 6 


0.0094 x 10° 


-8.1700 x 10° 
—8.1599 x 10* 

0.0101 x 10 6 


-8.1700 x 10° 
-8.1587 x 10 6 


0.0113 x 10° 


-8.1700 x 10° 
-8.1552 x 10* 

0.0148 x 10 6 


-8.1700 X 10° 
—8.1553 x 10* 


0.0147 x 10° 


— 8.1700 x 10° 
—8.1561 x 10* 

0.0109 x 10* 


-8.1700 x 10° 
-8.1555x10* 


0.0115 x 10° 


-8.1700 x 10° 
-8.1579x10* 

0A121 x 10* 


-8.1700x10* 
-8.1899 x 10* 


0*110 » 10 * 
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TABLE C-ll 

* TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 

HIGH ENERGY EQUATION SET 

f " MOON ASCENDING; 75 HOUR TRIP TIME 



3| 

'1 

>2 

F 

3 

H 

^NRT 

GEOCENTRIC 

ENERGY 

Pitchid Conic 

0.171S7 x 10* 

149.135 

12.077 

0 

48.0000 

-9.3a •'0 6 

Equation Sat 





-0.0788 

47.9316 

-9.293 •> 6 

Invars# 

0.17148 x 10* 

148.857 

11.720 




Error 

0.00019 x 10* 

0.278 

0.357 

0.0788 

0.0684 

0.0069 x 10 6 

Patchad Conic 

0.16261 x 10* 

89.674 

23.951 

0 

44.0000 

-9.3000 x 10 6 

Equation Sat 





-0.1531 

43.8808 

-9.2884 x 10 6 

Invars# 

0.16262 x 10° 

89.680 

22.987 




Error 

0.00001 x 10 s 

0.006 

0.964 

0.1531 

0.1192 

0.0116 x 10 6 

Patched Conic 

0.16164 x 10* 

43.413 

16.071 

0 

40.0000 

-9.3000 x 10* 

Equation Set 

-ft 




-0.1214 

39.8284 

-9.2809 x 10* 

Inverse 

0.16121 x 10® 

42.801 

15.160 




Error 

0.00043 x 10 s 

0.612 

0.911 

0.1214 

0.1716 

0.0191 x 10® 

Patched Conic 

0.16903 x 10* 

117.582 

13.807 

6.0000 

46.0000 

-9.3000 x 10 # 

Equation Set 





6.2199 

45.8856 

-9.2854 x 10* 

Inverse 

0.16807 x 10* 

116.776 

14.813 




Error 

0.00096 x 10 8 

0.806 


1.211 

0.2199 

0.1144 

0.0146 x 10® 

rasenea wontc 

0.16602 x 10* 

72.251 

12.286 

6.0000 

42.0000 

-9.3000 x 10* 

Equation Set 

ft 




6.2035 

41.8462 

-9.2809 x 10* 

Inverse 

0.16390 x 1C 8 

71.140 

13.712 




Error 

0.00112 x 10 s 

1.111 


1.427 

0.2035 

0.1638 

0.0191 x 10® 

Patched Conic 

0.12836 x 10* 

123.623 


12.0000 

52.0000 

-9.3000 x 10* 

Equation Set 





12.1275 

51.8369 

-9.2888 x 10* 

Inverse 

No Inverse 


■ 

■ 




Error 

Solution 


1 


0.1275 

0.0631 

0.0112 x 10® 

Patched Conic 

0.14110 x 10* 


60.518 

12.0000 

48.0000 

-9.3000 x 10® 

Equation Sat 





12.1458 

47.9188 

-9.2988x10* 

Inverse 

0.13884 x 10 8 

1 

■ 




Error 

0.00226 x 10* 

2.968 


0.1458 

0.0812 

0.0132 x 10® 

Patched Conic 

0.96894 x 10 7 

76.231 

154.880 

0 

54.0000 

-9.3000 x 10* 

Equation Set 





0.1121 

54.0840 

-9.3074 x 10* 

Inverse 

0.96663 x 10 7 

78.429 

156.712 




Error 

0.00041 x 10 7 

0.198 

0.852 

0.1121 

0.0840 

0.0074 x 10® 

Patched Conic 

0.97213 x 10 7 

53.522 

143.231 

2.0000 

52.0000 

-9.3000 x 10* 

Equation Set 





2.0841 

52.0651 

-0.8064 x 10* 

Inverse 

0.97028 x10 7 

56.283 

141.953 




Error 

a00416x 10 7 

1.741 


1.268 

0.0841 

0.0851 

0.0064 x 10® 


AVG 

0.1386 

0.1023 

i 

0.0121 x 10* 
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TABLE C-lll 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION 

HIGH ENERGY EQUATION SET 
MOON AT APOGEE; 60 HOUR TRIP TIME 


TO ANGULAR MOMENTUM 



Patched Conic 

0.85089 x 10 7 

123.081 

Equation Set 



Inverse 

0.85068 x 10 7 

122.941 

Error 

0.00021 x 10 7 

0.140 

Patched Conic 

0.84276 x 10 7 

164.371 

Equation Set 



Inverse 

0.84341 x 10 7 

166.079 

Error 

0.00065 x 10 7 

1.708 

Patched Conic 

0.86317 x 10 7 

132.241 

Equation Set 



Inverse 

0.862S9 x 10 7 

131.823 

Error 

0.00058 x 10 7 

0.418 

Patched Conic 

0.85657 x 10 7 

95.620 

Equation Set 



Inverse 

0.85605 x 10 7 

95.320 

Error 

0.00052 x 10 7 

0.300 

Patched Conic 

0.75480 x 10 7 

167.778 

Equation Set 



Inverse 

0.75495 x 10 7 

168.735 

Error 

0.00015 x 10 7 

0.957 

Patched Conic 

0.79979 x 10 7 

173.184 

Equation Set 



Inverse 

0.79953 x 10 7 

172.560 

Error 

0.00026 x 10 7 

0.624 

Patched Conic 

0.83314 x 10 7 

148.560 

Equation Set 



Inverse 

0.83223 x 10 7 

147.639 

Error 

0.00091 x 10 7 

0.921 

Patched Conic 

0.84474 x 10 7 

118.184 

Equation Set 

7 


Inverse 

0.84363 x 10 7 

117.480 

Error 

0.00111 x 10 7 

0.704 



39.070 

0.095 


^NRT 

^NRT 

GEOCENTRIC 

ENERGY 

0 

30.0000 

-2.8000 x 10 6 

—0.0210 

29.9809 

-2.7961 x 10 6 

0.0210 

0.0191 

0.0039 x 10 6 

3.0000 

32.0000 

-2.8000 x 10 6 

3.0216 

31.9763 

-2.7948 x 10 6 

0.0216 

0.0237 

0.0052 x 10 6 

3.0000 

30.0000 

-2.8000 x 10* 

3.0388 

29.9753 

-2.7947 x 10 6 

0.0388 

0.0247 

0.0053 x 10 6 

3.0000 

28.0000 

-2.8000 x 10 6 

0.0333 

27.9735 

-2.7944 x 10 6 

0.0333 

0.0265 

0.0056 x 10 6 

6.0000 

36.0000 

-2.8000 x 10 6 

6.0558 

35.9722 

-2.7930 x 10 6 

0.0558 

0.0267 

0.0070 x 10 6 

6.0000 

34.0000 

-2.8000 x 10 6 

6.0472 

33.9654 

-2.7917 x 10 6 

0.0472 

0.0346 

0.0083 x 10 6 

6.0000 

3Z0000 

-2.9000 x 10 6 

6.0501 

31.9643 

-2.7918 x 10* 

0.0501 

0.0357 

0.0082 x 10* 

6.0000 

30.0000 

-2.8000 x 10* 

6.0505 

29.9680 

-2.7927 x 10* 

0.0505 

0.0320 

0.0073 x 10* 

0.0398 

0.0279 

0.0064 x 10* 
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TABLE C-IV 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 

HIGH ENERGY EQUATION SET 
MOON AT APOGEE; 80 HOUR TRIP TIME 



ai 

i| 

■2 

^is,RT 

^NRT 

GEOCENTRIC 

ENERGY 

Patched Conic 

0.22120 x 10 s 

109.740 

44.074 

10.0000 

44.0000 

-8.1700 x 10 6 

Equation Sat 





10.0434 

43.5876 

-8.1270 x 10 6 

Invaraa 

0.21664 x 10 8 

107.039 

44.821 




Error 

0.00456 x 10 s 

2.701 

0.757 

0.0434 

0.4124 

0.0520 x 10 6 

Patchad Conic 

0.21537 x 10 s 

69.582 

44.075 

10.0000 

40.0000 

-8.1700 x 10 6 

Equation Sat 





10.0023 

39.5106 

—8.1190 x 10* 

Invaraa 

0.21233 x 10 s 

67.714 

44.073 




Error 

0.00304 x 10 s 

0.232 

0.003 

0.0023 

0.4894 

0.0610 x 10 6 

Patchad Conic 

0.13803 x 10 8 

51.911 

86.116 

14.0000 | 

48.0000 

-8.1700 x 10 6 

Equation Sat 





14.3575 

48.1806 

-8.1708 x 10 6 

Invaraa 

0.13225 x 10* 

40.670 

90.363 




Error 

0.00578 x 10 s 

11.241 

4.247 

0.3575 

0.1806 

0.0008 x 10* 

Patchad Conic 

0.18616 x 10 s 

32.749 

56.441 

14.0000 

40.0000 

-8.1700 x 10 6 

Equation Sat 





14.2389 ! 

39.7658 

-8.1367 x 10 6 

Invaraa 

0.18374 x 10 8 

29.813 

56.907 




Error 

0.00242 x 10* 

0.064 

0.466 

0.2389 

0.2342 

0.0333 x 10 6 

Patchad Conic 

0.12871 x 10* 

101.125 

95.887 

20.0000 

54.0000 

-8,1700 x 10 6 

Equation Sat 





20.4039 

54.2541 

-8.1732 x 10 6 

Invaraa 

No Invaraa 







Error 

Solution 




0.4039 

0.2541 

0.0032 x 10 6 

Patchad Conic 

0.16331 x 10 s 

48.160 

69.479 

20.0000 

44.0000 

-8.1700 x 10 6 

Equation Sat 





20.4081 

44.0395 

-8.1555 x 10 6 

Invaraa 

0.16616 x 10 8 

59.160 

66.332 




Error 

0 00585 x10 8 

11.000 

1 3.147 

1 

0.4081 

0.0395 

0.0146 x 10 6 


AVG j 

0.2424 

0.2684 

0.0258 x 10® 
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TABLE C-V 


t 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 

HIGH ENERGY EQUATION SET 
MOON AT APOGEE; 80 HOUR TRIP TIME 



at 

*1 

*2 

Asirt 

^NRT 

GEOCENTRIC 

ENERGY 

Patched Conic 

0.15186 x 10 S 

154.571 

87.655 

0 

-72.0000 

-8.1700 x 10 6 

Equation Set 




0.9825 

71.5678 

-8.0956 x 10 6 

Inverse 

0.15416 x 10 8 

157.413 

82.853 




Error 

0.00230 x 10 8 

0.158 

4.802 

0.9825 

0.4322 

0.0744 x 10® 

Patched Conic 

0.21048 x 10 s 

169.435 

22.472 

0 

69.0000 

-8.1700 x 10® 

Equation Set 




0.25)89 

68.6042 

-8.1287 x 10® 

Inverse 

0.20610 x 10 8 

167.212 





Error 

0.00438 x 10 8 

2.223 


0.2989 

0.2958 

0.0S13 x 10® 

Patched Conic 

0.21478 x 10 S 

172.181 

14.413 

5.0000 

69.0000 

-8.1700 x 10® 

Equation Set 




5.4090 

68.5826 

-8,1240 x 10® 

Inverse 

0.21209x10® 

171.242 

14.916 




Error 

0.00269 x 10® 

0.939 

0.503 

0.4090 

0.4174 

0.0450 x 10® 

Patched Conic 

0.19950 x 10® 

150.336 

33.350 

5.0000 

66.5000 

-8.1700 x 10® 

Equation Set 




4.8211 

66.1798 

-8.1392 x 10® 

Inverse 

0.19266 x 10® 

148.386 

39.409 




Error 

0.00684 x 10® 

1.950 

6.059 

0.1789 

0.3202 

0.0338 x 10® 

Patched Conic 

0.17941 x 10® 

169.343 

59.240 

10.0000 

72.0000 

-&1700 x 10® 

Equation Set 




10.7819 

70.9716 

-8.0068 x 10* 

Inverse 

0.17018x10® 

167.908 

62.856 




Error 

0.00923 x 10® 

1.435 

3.616 

0.7819 

1.0284 

0.1632 x 10* 

Patched Conic 

0.16050 x 10® 

116.810 

66.759 

10.0000 

58.4000 

-8.1700 x 10® 

Equation Set 




9.7810 

58.1513 

-8.1364 x 10® 

Inverse 

No Inverse 






Error 




0.2190 

0.2487 

0.0308 x 10® 

Patched Conic 

0.20726 x 10® 

144.973 

19.262 

15.0000 

65.0000 

-8.1700 x 10® 

Equation Set 




15.8146 

64.2683 

-8.0784 x 10® 

Inverse 

0.20046 x 10® 

143.942 

23.235 




Error 

0.00680 x 10® 

1.031 

3.973 

0.8146 

0.7317 

0.0916 x 10® 

Patched Conic 

0.14693 x 10® 

143.782 

92.059 

20.0000 

71.0000 

-8.1700 x 10® 

Equation Set 




20.6139 

69.6567 

-8.0008 x 10* 

Inverse 

0,13539 x 10 s 

140.093 

102.1224 




Error 

0.01164 x 10® 

3.689 

10.063 

0.6139 

1.3443 

0.1692 x 10® 

Patched Conic 

0.20047 x 10® 

121.511 

13.585 

20.0000 

57.0000 

-8.1700 x 10* 

Equation Set 




20.8720 

58.2725 

-8.0774 x 10® 

Inverse 

0.19404 x 10® 

120.633 

18.519 




Error 

0.00643 x 10® 

0.878 

4.934 

0.8720 

0.7274 

0.0926 x 10* 

Patched Conic 

0.11248x10® 

97.537 

140.941 

26.0000 

57.0000 

—8.1700 x 10* 

Equation Set 




24.4829 

58.5375 

-8.1286 x 10* 

Inverse 

0.11046x10* 

96.542 

145.441 




Error 

0.00202 x 10® 

0.995 

4.600 

0.5171 

0.1625 

0.0414 x 10® 

AVQ 

0.5888 

0.6709 

0.0793 x 10* 


AVQ 
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TABLE C-VI 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 

HIGH ENERGY EQUATION SET 
MOON AT APOGEE; 80 HOUR TRIP TIME 


a i 




RT 


V N RT 


GEOCENTRIC 

ENERGY 


Patched Conic 
Equation Set 
Inverse 
Error 


0.21501 x 10° 

0.22247 x 10 8 
0.00746 


111.348 

112.275 

0.927 


51.354 

48.533 

2.821 


4.0000 

4.3792 

0.3792 


54.0000 

53.8165 

0.1835 


-8 1700 x 10 6 
-8.1683 x 10 € 


0.0017 x 10' 


Patched Conic 
Equation Set 
Inverse 
Error 


0.18270 x 10* 


8 


0.18202 x 10 
0.00068 x 10 


93.661 

93.539 

0.122 


61.484 


60.675 

0.809 


4.0000 

4.1258 

0.1258 


48.0000 

47.3723 

0.6277 


-8.1700 x 10° 
-8.1141 x 10* 


0.0559 x 10' 


Patched Conic 
Equation Set 
Inverse 
Error 


0.22034 x 10° 
0.23301 x 10 8 


0.01267 x 10 J 


8 


100.692 

102.050 

1.358 


46.322 

41.134 

5.188 


6.0000 

7.0660 

1.0660 


48.0000 

47.5005 

0.4995 


-8.1700 x 10° 
-8.1568 x 10* 

0.0132 x 10 6 


Patched Conic 
Equation Set 
Inverse 
Error 


0.21039 x 10 8 

0.21706 x 10* 
0.00667 x 10 s 


91.465 

92.306 

0.841 


47.984 

43.882 

4.102 


6.0000 

6.8737 

0.8737 


44.0000 

43.1564 

0.8436 


-8.1700 x 10 6 
-8.1134 x 10 6 

0.0566 x 10* 


Patched Conic 
Equation Set 
Inverse 
Error 


0.85754 x 10 


0.86579 x 10 ' 
0.00825 x 10 7 


57.918 

58.805 

0.887 


159.124 

154.265 

4.859 


12.0000 

11.0533 

0.9467 


48.0000 

47.9265 

0.0735 


-8.1700 x 10 6 
-8.1397 x 10 6 

0.0303 x 10 6 


Patched Conic 
Equation Set 
Inverse 
Error 


0.81972 x 10 ' 


0.82075 x 10 7 
0.00103 x 10 7 


28.021 

28.284 

0.263 


167.410 

169.273 

1.863 


12.0000 

11.6782 

0.3218 


40.0000 

40.2808 

0.2808 


-8.1700 x 10 6 
-8.2026 x 10 6 

0.0326 x 10 6 


Patched Conk 
Equation Set 


Error 


0.87779 x 10' 


0.87744 x 10 7 
0.00035 x 10 7 


74.928 

74.893 

0.035 


160.089 

163.885 

3.796 


24.0000 

23.1392 

0.8608 


54.0000 

54.5368 

0.5368 


-8.1700 x 10° 
-8.2361 x 10* 

0.0661 x 10 6 


I 

l 


Patched Conk 
Equation Set 


Error 


Patched Conk 
Equation Set 


Error 


0.93289 x 10 


0.93618 x 10 7 
0.00329 x 10 7 


0.12203 x 10° 


0.12118 x 10° 
0.00085 x 10 8 


43.986 

44.477 

0.491 


131.307 

133.776 

2.469 


24.0000 

23.7883 

0.2117 


88.743 

88.441 

0.302 


101.872 

102.832 

0.960 


36.0000 

35.8902 

0.1098 


40.0000 

41.1382 

1.1382 


54.0000 

54.0864 

0.0864 


-8.1700 x 10 6 
-8.3020 x 10 6 

0.1320 x 10* 


-8.1700 x 10*i 
-8.1316x10* 


0.0116 x 10® 


Patched Conk 
Equation Set 


0.14466 x 10° 


0.14773 x10 8 
0.00307 x 10* 


69.777 

70.734 

0.967 


79.386 

78.974 

0.412 

H 


36.0000 

36.3200 

0.3200 


40.0000 

40.8693 

0.8693 


-8.1700 x 10* 
-8.2639 x 10* 

0.0839x10* 


Error 


AVG 


0.5216 


0.8109 


(MMMx id* 




































































TABLE C-VII 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 

GENERAL EQUATION SET 
MOON ASCENDING; 75 HOUR TRIP TIME 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 


Patched Conic 
Equation Set 
Inverse 
Error 




0.17167 k 10® 

0.17116 x 10® 
0.00051 x 10® 



0.16164 x 10 s 


0.16110 x 10° 
0.00054 x 10® 


0.16903 x 10® 

0.16808 x 10® 
0.00095 x 10® 



116.780 

0.802 


0.12836 x 10” 

0.12416 x 10® 
0.00420 x 10* 


0.14110x 10* 

0.13886 x 10* 
0.00224 x 10* 


0.096694 x 10* 

0.096636 x 10* 
0.000068x10* 


0.097213 x 10 s 

0.097960 x 10* 
0.000737 x 10* 


12.0000 

12.1950 


40.0000 

39.8839 


46.0000 

45.8924 


52.0000 

51.9827 


48.0000 

47.9186 


GEOCENTRIC 

ENERGY 


-9.3000 x 10 6 
-9.2875 

0.0125 x 10® 


Patched Conic 

0.16261 x 10® 

89.674 

23.951 

0 

44.0000 

-9.3000 x 10® 

Equation Set 

0.16088 x 10® 



0.3429 

43.8836 

-9.2826 x 10® 

Inverse 

88.095 

■ 




Error 

0.00173x10® 

1 

1.579 

I 

0.3429 

0.1164 

0.0174 x 10® 


-9.3000 x 10° 
-9.2870 x 10® 


0.0130 x 10° 


-9.3000 x 10° 
-9.2860 x 10® 


0.0140 x 10° 


Patched Conic 

0.16502 x 10® 

72.251 

12.285 

6.0000 

42.0000 

-9.3000 x 10® 

Equation Set 




6.1133 

42.5215 

-9.3970 x 10® 

Inverse 

0.16840 x 10* 

71.681 

12.725 




Error 

0.00338x10® 

0.570 

0.440 

0.1133 

0.5215 

0.0970 x 10® 


-9.3000 x 10° 
-9.2914 x 10® 

0.0086 x 10® 


-9.3000 x 10® 
-9.2888 x 10® 

0.0132 x 10® 



-9.3000 x 10” 
-9.3076 x 10* 

0.0076 x 10* 


-9.3000 x 10* 
-9.3094 x 10* 



































































TABLE C-VII I 


TRANSLUNAR TRAJECTORIES USING THE PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 


GENERAL EQUATION SET 
MOON AT PERIGEE; 90 HOUR TRIP TIME 





»2 

^NRT 

x 

GEOCENTRIC 


a l 


a nrt 

ENERGY 

Patched Conic 

0.26270 x 10 8 

168.446 

4.561 

0 

74.5000 

-11.4700 x 10 6 

Equation Sat 





0.1687 

73.845? 

-1 1.4500 x 10 6 

Inverse 

0.26123 x 10 s 

171.190 

4.252 




Error 

0.00147 x 10 8 

2.744 

0.309 

0.1687 

0.6541 

0.0200 x 10 6 

Patched Conic 

0.23825 x 10 s 

91.765 

22.871 

0 

65.0000 

-11.4700 x 10 6 

Equation Set 





0.1494 

64.5692 

—11.4523 x 10* 

Inverse 

0.23616 x 10 s 

91.428 

23.620 




Error 

0.00209 x 10 s 

0.337 

0.749 

0.1494 

0.4308 

0.0177 x 10* 

Patched Conic 

0.23707 x 10 s 

18.169 

6.801 

0 

56.0000 

-11.4700 x 10 6 

Equation Set 





0.1291 

55.6671 

—11.4510 x 10* 

Inverse 

0.22743 x' 8 

17.619 

7.199 




Error 

0.00964 x 1u“ 

0.550 

0.398 

0.1291 

0.3329 

0.0190 x 10 6 

Patched Conic 

0.23848x 10 s 

168.369 

31.391 

7.0000 

75.0000 

-11.4700 x 10 6 

Equation Set 





7.2880 

74.5898 

-11.4502 x 10* 

Inverse 

0.22743 x 10 8 

172.645 

miss 




Error 

0.01105 x10 s 

3.286 

6.764 

0.2880 

0.4102 

0.0198 x 10* 

Patched Conic 

0.25000 x 10 s 

98.910 

10.398 

7.0000 

66.0000 

-11.4700 x 10 6 

Equation Set 





7.5783 

65.4860 

-11.4421 x 10* 

Inverse 

0.24677 x 10 s 

99.062 

13.242 




Error 

0.00328 x10 s 

0.848 

2.844 

0.5783 

0.5140 

0.0279 x 10* 

Patched Conic 

0.22462x10* 

23.316 

25.585 

7.0000 

57.0000 

—11.4700 k 10* 

Equation Set 





7.1265 

56.7654 

-11.4543x 10* 

Inverse 

0.22243 x 10* 

22.590 

26.700 




Error 

0.00219 x 10* 

0.726 


1.115 

0.1255 

0.2346 

0.0157 x 10* 

Patched Conk 

0.17428x10* 

119.455 

71.605 

14.0000 

72.0000 

-11.4700x10* 

Equation Set 





14.2612 

71.8253 

-11.4526 x 10* 

Inverse 

Error 

No Inverse 




0.2812 

0.1747 

0.0175 x 10* 

Patched Conk 

0.19028 x 10* 

96.847 

58.614 

14.0000 

eaoooo 

-11.4700 x 10* 

Equation Set 





13.9309 

87.7128 

-11.4592x10* 

Inverse 

0.19099x10* 

96.920 

67.802 




Error 

0.00071 x 10* 

0.073 

0.812 

0.0691 

0.2872 

0.0109 x 10* 

Patched Conk 

0.17658 x 10* 

19 

65.996 

14.0000 

64.5000 

-11.4700x10* 

Equation Set 
Inverse 

at — lauiaflae 

no invifis 

fjitf 



14.0769 

64.4401 

-11.4642x10* 

Error 


1 



0.0769 

0.0699 

0.0068x10* 

Patched Conk 

0.10844x10* 

92.149 

160.000 

6.7600 

72.3670 

-11.4700x10* 

Equation Set 





■ 

72.6071 

-11.4661 x 10* 

Inverse 

0.10092x10* 

93.163 





Error 

0.00039 x 10* 

1.014 

0.462 


0.4201 

0.0181 x 10* 

Patched Conk 

aioon x io* 

127.429 

190.000 

6.0000 


-114700x10* 

Equation Set 





8.7679 

■ i 

-114077x10* 

Inverse 

0.10993x10* 

127.712 

191.670 


wKHm 


Error 

0.00010 x 10* 

0.293 


1.670 

0.2321 


04277 x 10* 

Patched Conk 

0.10617x10* 

16.799 

100.000 

4.0000 

66.6620 

-114700x10* 

Equation Set 





3.9026 

66.8796 

-114867x10* 

Inverse 

0.19623 x 10* 

19.990 

161.202 



0.0167x10* 

Error 

OOOOOOx 10* 

1.134 


1.202 

0.0974 

04276 



AVG 

0.1977 

04724 

0410x10* 
























































TABLE C-IX 



SUMMARY OF ERRORS FOR TRANSLUNAR TRAJECTORIES USING THE 
PARABOLIC APPROXIMATION TO ANGULAR MOMENTUM 


INITIAL 

CONDITIONS 

NO. 

OF 

RUNS 

EQN. 

SET* 

ERRORS 


LAT. 

LONG. 

I GEO. ENERGY 

DEG. 

DLG. 


H pc - 100 N. Ml. 

Hpg - 60 N. Ml. 

E e - -8.17 x 10 6 -FT 2 /SEC 2 
MOON AT APOGEE 

12 

H 

AVG. 

0.1213 

0.0900 

0.0110 X 10 6 

MAX. 

■ 

■ 

m 

H pG = 100 N. Ml. 

Hpg- 60 N. Ml. 

Eg - -9.30 x 10 6 FT 2 /SEC 2 
MOON AT MID DISTANCE 

9 

H 

AVG. 

■ 

■ 

■i 

MAX. 

■ 

■ 

B 

HpG - 100 N. Ml. 

Hpg - 60 N. Ml. 

E e - -2.8 x 10 s FT a /SEC 2 
MOON AT APOGEE 

8 

H 

AVG. 

■ 

■ 

0.0064 x 10 6 

MAX. 

0.0558 

0.0357 

0.0083 x 10 6 

HpQ - 10,000 N. Ml. 
Hpg- 60 N. Ml. 

6 

H 

AVG. 

0.2424 

0.2884 

0.0258 x 10 6 

E e - —8.17 x 10* FT 2 /SEC 2 
MOON AT APOGEE 

MAX. 

0.4081 

0.4894 


Hpg- 100 N. Ml. 

Hpg- 10,000 N. Ml. 

E e - -8.17x10* FT 2 /SEC 2 
MOON AT APOGEE 

10 

H 

AVG. 

H 

0.5682 

0.0700 x 10 6 

MAX. 

0.9825 

1.3443 

B1 

Hpg - 10,000 N. Ml. 

Hpg - 10,000 N. Ml. 

E e - -8.17 x 10 6 FT 2 /SEC 2 
MOON AT APOGEE 

10 

H 

AVG. 

0.5216 

0.5109 


MAX. 

1.0660 

1.1382 

0.1320 x 10 6 

Hpg - 100 N. Ml. 

Hpg - 60 N. Ml. 

E e - -9.30 x 10 6 FT 2 /SEC 2 
MOON AT MID-DISTANCE 

9 

G 

AVG. 

0.1581 

0.1422 

0.0193 xIO 6 

MAX. 

0.3429 

0.5417 

0.0970 x10 s 

Hpg - 100 N. ML 
Hpg -SON. Ml. 

E e - -11.47 x 10* FTVSEC 2 
MOON AT PERIGEE 


G 

AVG. 

0.1877 

0.3724 

0.0183 x10 s 

MAX. 


g 

0.0279 x10 s 


*H REFERS TO THE HIGH ENERGY EQUATION SET 
G REFERS TO THE GENERAL EQUATION SET 

HpQ -PERIGEE HEIGHT 
Hpg -PERISELENE HEIGHT 
E e - GEOCENTRIC ENERGY 


















































































